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neutrino oscillations taking account of the matter effects in the earth. With 



i-S^ \ the hierarchy among the vacuum mass eigenvalues /i^ ^ /U^ > /uf , the solution 

Q-i! 

Q^' of the atmospheric neutrino problem depends on = fi'^ — fil and the 13 and 



23 mixing angles (p and ip. Whereas the sub-GeV atmospheric neutrino data 
imposes only a lower limit on 631 > 10^^ eV'^, the zenith angle dependent 



X 

I suppression observed in the multi-GeV data limits ^31 from above also. The 

allowed regions of the parameter space are strongly constrained by the multi- 
GeV data. Combined with our earlier solution to the solar neutrino problem 
which depends on 621 = /^^ — /^i and the 12 and 13 mixing angles u> and (/), 
we have obtained the ranges of values of the five neutrino parameters which 
solve both the solar and the atmospheric neutrino problems simultaneously. 
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I. INTRODUCTION 



In recent years, deep underground detectors have been measuring the flux of neutrinos 
produced in the atmosphere. (Unless exphcitly stated, we call neutrinos and anti-neutrinos 
collectively as neutrinos). These neutrinos are produced from the decays of vr^ and 
mesons which, in turn, are produced by cosmic ray interactions with the atmosphere. One 
expects the number of muon neutrinos to be twice the number of electron neutrinos. The 
predictions of detailed Monte Carlo calculations conflrm that the ratio of the flux of muon 
neutrinos $1,^ to the flux of electron neutrinos (^^^ is about 2 [|T],3. The predicted neutrino 
fluxes from different calculations differ signiflcantly from each other (by as much as 30%), 
but the predictions for the ratio of the fluxes ^^^/(^y^ are in good agreement with each 
other. The large water Cerenkov detectors Kamiokande and 1MB have measured this 
ratio and have found it to be about half of what is predicted [^. The experimental results 
are presented in the form of a double ratio 



Kamiokande collaboration have presented their results for neutrinos with energy less 



(multi-GeV data) §. For the sub-GeV data, R = OMto fa ± 0.05 and for the multi-GeV 
data, after averaging over the zenith angle, R = 0.57lo;o7 ± 0.07. The value of R has no 
signiflcant zenith angle dependence for the sub-GeV data. However, for the multi-GeV data, 
R is smaller for large values of zenith angle (upward going neutrinos) and is larger for small 
values of zenith angle (downward going neutrinos). Kamiokande have analyzed their data 
assuming that the smaller observed value of R is caused by neutrino oscillations. They have 
done two independent analyses, one assuming two flavor oscillations between ^ z/g and 
the other assuming two flavor oscillations between ^ Ur- For both the cases, they obtain 
a mass-squared difference 6171^ ~ 10~'^eV^ and a mixing angle nearly 45°. Since the upward 
going neutrinos travel large distances inside the earth before entering the detector, matter 



R 




(1) 



than 1.33 GeV (sub-GeV data) |^ and for neutrinos with energy greater then 1.33 GeV 
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effects may be important for these, especially at higher energies Therefore one must 

take matter effects into consideration while analyzing <-> z/g oscillations. 

Atmospheric neutrino problem was analyzed in the three flavor neutrino oscillation frame- 
work previously. In Ref. 0] the sub-GeV data were analyzed under the assumption that one 
of the mass differences is much smaller than the other. The matter effects due to the passage 
of neutrinos through the earth were included and the allowed values of neutrino parameters 
were obtained. Various other authors have analyzed, in the context of three flavor oscilla- 



tions, accelerator and reactor data in conjunction with the sub-GeV atmospheric data ||1C 



or the multi-GeV data with zenith angle dependence included [|Tl|]. However, in both these 
cases the earth matter effects were not taken into account. 

Several authors have attempted a simultaneous solution of the solar and atmospheric 
neutrino problems in three flavor oscillation scenarios. The solution in Ref. fl^ 



assumes 



maximal mixing between all the three flavors and derives the constraints on the mass differ- 
ences. The solution thus obtained restricts the mass difference to be very small (~ 10^^° eV"^) 
and is somewhat fine-tuned. Other solutions have assumed the mass hierarchy that was con- 
sidered in Ref. and obtained the allowed regions in the neutrino parameter space [T^-|TB 



However, all these analyses were based on the sub-GeV atmospheric neutrino data. The 
multi-GeV data and its zenith angle dependence were not included. In a recent analysis of 
the multi-GeV data earth matter effects and the zenith angle dependence were taken into 



account [T^. However, this analysis considered the atmospheirc neutrino problem only. 

In this paper, we analyze the Kamiokande multi-GeV data in the framework of three 
flavor neutrino oscillations. The solar neutrino data and the sub-GeV data of Kamiokande 



were analyzed earlier in the same framework |jT6[. We only assume that, of the two mass 
differences in the problem, the smaller one is chosen to solve the solar neutrino problem and 
the larger one will be relevant for the atmospheric neutrino problem. We do not, apriori, 
make any assumption about the form of the mixing matrix. Our starting point is the same 
as that in Ref. . We include the matter effects due to the passage of neutrinos through the 
earth. We take the allowed values of neutrino parameters from our earlier work as inputs. 



The zenith angle dependent suppression of multi-GeV data places strong further constraints 
on the parameter space allowed by the solar neutrino problem and the sub-GeV atmospheric 
neutrino data. In section II, we define our theoretical framework and derive the oscillation 
probabilities. In section III, we describe the comparison of the theoretical predictions with 
the Kamiokande data and determine the allowed regions of the neutrino parameters. In 
section IV we discuss our results and present our conclusions. A brief discussion of the 



recent LSND result 18 is also included. 



II. THEORY 

Here we describe the three fiavor neutrino mixing and calculate the probability for a 
neutrino produced as z/^ to be detected as up, where a and /5 are fiavor indices. The fiavor 
eigenstates are related to the mass eigenstates by 











= 








/3_ 



where the superscript 'v' denotes vacuum. The unitary matrix V" can be parametrized as 

U"" = U^\iIj) X UP''''"' X U^\(j)) X U^\uj), (3) 

where W^{6ij) is the two fiavor mixing matrix between the ith and jth mass eignestates with 
the mixing angle 9ij. For simplicity, we neglect the CP violation and set JJ'p'^"-^'' = J. 
In the mass basis the mass-squared matrix is diagonal and can be taken to be 
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^31 



where 621 = ^-^^ ^31 = fif. Without loss of generality we can take /is > yU2 > yUi 

so that both and ^31 are positive. For extreme relativistic neutrinos, the oscillation 
probability depends only on the mass-squared differences so we ignore the first term ^fl. 
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In the flavor basis, the mass- squared matrix is non-diagonal and is given by 

Ml = U^Ml iU'')^ . (5) 

For the propagation of neutrinos through the earth, we need to take the matter effects into 
account. The charged current scattering between electrons and z/g induces an effective mass- 
squared term for i/g which is of the form A = 2\^GFNeE, where Ne is the number density 



of electrons and E is the neutrino energy |TP[ . This term is present only for the e — e element 
in the flavor basis so the mass-squared matrix including the matter effects is 



A 





(6) 











= f/™ 


m 






m 
^3 



is a hermitian matirx and can be diagonalized by unitary matrix t/*", which relates the 
flavor eigenstates to matter dependent mass eigenstates 



(7) 



where the superscript 'm' stands for matter. We denote the matter dependent mass eigen- 
values as mi, 1712 and 7713. The matter dependent mixing matrix f/"* can be parametrized in 
terms of three mixing angles in a manner similar to that of V" as given in eq. (|^) . The mat- 
ter dependent mass eigenvalues and mixing angles can be obtained, in terms of the vacuum 
parameters and A, by solving the eigenvalue problem of in eq.(^ 

The distance scales and the energy scales in the solar neutrino problem and the atmo- 
spheric neutrino problem are very different. Therefore, one needs two distinct mass scales 
to solve the solar and the atmospheric neutrino problems simultaneously. To satisfy the 
constraints from all the three solar neutrino experiments, one must choose ^21 ~ 10~^ eV'^, 
which is roughly the matter term A for the solar neutrinos due to their passage through the 
sun. As we will shortly see, the sub-GeV atmospheric neutrino data impose a constraint 



^31 > 10~^ eV'^. Hence S^i ^ 621- If ^21 ~ 10~^ eV'^, the oscillation length corresponding 
to it, even for the minimum of the atmospheric neutrino energies, is of the order of the di- 
ameter of the earth. In the expression for the oscillation probability, 621 can be set to zero. 
Therefore the oscillations are dependent on only one mass difference ^31 in the atmospheric 
neutrino problem. 

In the approximation of neglecting ^21, it is straightforward to show that the oscillation 
probability is independent of the mixing angle u. Including matter effects does not alter this 
conclusion. Note, however, that neglecting 621 does not reduce the problem to an effective 
two flavor mixing. The three flavor nature of the problem is reflected by the fact that the 
oscillation probability is a function of the mass difference ^31 and two mixing angles and 
ip. In the case of an effective two flavor mixing, the oscillation probability is dependent on 
one mixing angle only. When the matter effects are included the mixing angle ip remains 
unaffected but the angle becomes matter dependent 0, 

^31 sin 20 

tan20^ = - — -. (8) 

631 cos 20 — ^4 

The mass eigenvalue of remains (actually it is of the order of ^21 which we are neglecting 
here). The other two matter dependent mass eigenvalues are given by 

o 1 



rrir, = - 

2 2 

2 1 

^ 2 



(^31 + A)- V(53i cos 20 - AY + (531 sin 20)^ 



(^31 + A) + V (^31 cos 2<p-Af + (^31 sin 20)' 



(9) 
(10) 



Equations (^, (P) and (|TOD are valid for neutrinos. For anti- neutrinos, we get a similar set 
of formulae with A replaced by —A. 

The neutrinos produced in the atmosphere enter the earth after travelling through the 
atmosphere for about 20 Km and finally reach the detector after travelling through the 
earth. The distance travelled through the earth is a function of the zenith angle. For the 
five bins considered by Kamiokande , the average values of the cosine of the zenith angle 
are —0.8, —0.4, 0.0, 0.4, 0.8 and the average distances travelled through the earth are 10210, 
5137, 832, 34, 6 Km respectively H . 



A neutrino of flavor a, produced in the atmospliere at time t = 0, propagates through 
the atmosphere as a hnear combination of the vacuum mass eigenstates. If the neutrino 
enters earth at time t = ti, its state vector at that time can be written as 

\^M) = Y.U:.exp|^-^^yu^). (11) 

Reexpressing the vacuum mass eigenstates in terms of flavor states, we have 

l^aih)) =Y:K,exp Ef^A7kA). (12) 

After entering the earth, the neutrino propagates as a linear combination of the matter 
dependent mass eigenstates. We take the earth to be a slab of constant density. At the time 
of detection t = t^, the state vector takes the form 

Hence the amplitude for the neutrino produced as flavor a at t = to be detected as a 
neutrino of flavor 13 at time td is given by 

{M^M) = T.T.T.U:.U::UT,U^*exp (^-.^) exp (-^^^^^^^^) ■ (14) 

The probability of oscillation Pap is given by the modulus square of the above amplitude. 
If td — ti is set equal to zero (that is if the total time of travel is equal to the time of 
travel through the atmosphere) then the expression in eq. ([I^ ) reduces to the simple vacuum 
oscillation amplitude. The same is true if the matter effects are ignored; i.e. if U"^ = V" 
and rrii = fii. 

For bins 1,2 and 3, the distance travelled in earth is much larger than the distance 
travelled in the atmosphere. Therefore ti is much smaller than td for these bins and can be 
neglected. Neglecting ti, simplifies the expressions for oscillation probability and we obtain 
the expressions derived earlier in Ref. 0]. However, for bins 4 and 5, the distance of travel in 
atmosphere is comparable to that in earth. Therefore ti is of the same of order of magnitude 
as td and can not be neglected. Keeping ti 7^ in eq.(|Tl|) properly takes into account the 
non-adiabaticity in the abrupt change in density when the neutrino enters earth. 
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III. CALCULATION AND RESULTS 
A. Sub-GeV Data 

First we briefly describe our analysis of the sub-GeV data and highlight the contrast 
between it and the analysis of the multi-GeV data. Matter effects are unimportant for the 
sub-GeV data. If the earth is taken to be a slab of density 5.5 gm/cm^, the matter term A 
for the sub-GeV neutrinos is less than 3.8 x 10"'^ eV"^. As we will shortly see, the sub-GeV 
data sets a lower limit on 631 > 10~^ eV"^. Hence the matter effects can be neglected and the 
expressions for Pa/B in the sub-GeV analysis are simply the vacuum oscillation probabilities 

P'a, = {u:,u;,y + {u:,u;,y + {u:,u;,y + 2 u:,u:,u;,u;, COS (2.53^) + 

2 Ul,Ul,Ul,Ul, cos (2.53^^ + 2 Ul,Ul,Ul,Ul, cos ^2.53^^ , (15) 

where d is the distance of travel in meters, (5's are the mass differences in eV"^ and E is the 
neutrino energy in MeV. Because we have neglected the CP violating phase, the oscillation 
probability for the anti-neutrinos is the same as that for the neutrinos. Since is very 
small, the cosine term containing it in eq.(|15D can be set equal to 1. The other two cosine 
terms are dependent on ^si (and ^32 — ^si), the neutrino energy and the distance of travel 
which is related to the zenith angle. As mentioned earlier, the double ratio R defined in 
eq.(|l]) does not have any zenith angle dependence for the sub-GeV data. One can account for 
this if it is possible to replace the distance dependent terms in eq. ( P^D by their average values. 
This replacement is possible only if the average distance travelled contains many oscillation 
lengths. The above condition sets a lower limit on the mass difference ^31 > 10~^ eV^ 
||T6| . Note that this lower limit is consistent with the approximation ^31 ^ ^21 ~ 10~^ eV'^, 
which was made so that both solar and atmospheric neutrino problems could be solved 
simultaneously. 

P^p become independent of 531 and the neutrino energy when the distance dependent 
terms are replaced by their average values. They are simply functions of the mixing angles 
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and if) alone. In this approximation, the expression for the double ratio R can be written 



as [ra 



pO J ]_po 

where tmc is the Monte Carlo expectation of the ratio of number /i-like events to the 
number of e-like events. From the sub-GeV data of Kamiokande, we find tmc = 1-912 and 
R = 0.60l[i:[|^±0.05 |,|1|]. We take the allowed values of from our analysis of solar neutrino 
data which restricts to be in the range < < 50° [Q. We find the allowed region in 



the (j) — 4' plane by requiring the theoretical value calculated from eq. (piq) to be within la 
or 1.6a uncertainty of the experimental value. In our previous analysis we restricted the 
range of -0 to be < -0 < 45°. However, given the assumptions we made about vacuum 



mass eigenvalues, the most general possibility is to allow ip to vary between and 90 |^ 
The new region allowed by the sub-GeV data, where was restricted by the solar neutrino 
data and ip is allowed to vary between to 90", is shown in Figure 1. The region between 
the solid lines is the parameter space which satisfies the experimental constraints at la level 
whereas the region between dashed lines satisfies the experimental constraints at 1.6a level. 
One important point to be noted in this analysis is that the sub-GeV data place only lower 
bound on ^31. The allowed region in — -0 plane is quite large. 

B. Multi-GeV data 



The multi-GeV data of Kamiokande have been presented for five zenith angle bins in 
Ref. 0. For each of these bins, the observed numbers of electron-like events and muon-like 
events and their Monte Carlo expectations (without neutrino oscillations) have been given. 
From these one can calculate two sets of ratios 

'■'-c = (§) (17) 

V^^e/ MC 

<>^s=[j^] , (18) 

V^^e / obs 



and the set of double ratios 



^obs — ''^^obs/^^MC (19) 

for each bin i = 1,2,...,5. We summarize the multi-GeV data of Kamiokande in Table I. 

The /i-like events are subdivided into fully contained (FC) and partially contained (PC) 
events whereas all the e-like events are fully contained. The efficiency of detection for each 
type of event is different and is a funtion of the neutrino energy. Thus we have three detection 
efficiencies epQ{E), ep(j{E) and s^{E). We obtained these efficiencies from Kamiokande 
Collaboration . 



The expected number of /x-like and e-like events, in the absence of neutrino oscillations, 
is given by 

NlUi.c. = J [€SE)aiE) + <l>l^iE)aiE)] (^^(E) + ^^^(E)) dE (20) 
I [<PliE)aiE) + <Pl{E)aiE)] e\E)dE (21) 



\M.C. 



where 0's are the fluxes of the atmospheric neutrinos at the location of Kamiokande. These 
are tabulated in Ref. [|I| as functions of the neutrino energy (from E = l.Q GeV to E = 100 
GeV) and the zenith angle, a and a are the charged current cross sections of neutrinos and 
anti-neutrinos respectively with nucleons. The cross section is the sum of the quasi-elastic 
scattering and the deep inelastic scattering (DIS). The values for quasi-elastic scattering 



are taken from Gaisser and O'Connel p2[ and those for the DIS are taken from Gargamelle 
data [^]. In calculating the DIS cross section, we took the lepton energy distribution to 
be given by the scaling formula (which is different for neutrinos and anti-neutrinos) and 
integrated da/dEi^p from the minimum lepton energy Emm = 1-33 GeV to the maximum 
lepton energy E^ax = Ey — m,r- The maximum lepton energy is chosen by deflning DIS to 
contain at least one pion in addition to the charged lepton and the baryon. The differences 
in the flducial volumes and exposure times for fully contained and partially contained events 



have been incorporated into the detection efficiency e^(j{E). From equations ( [2(]| ) and ( pi]) 
we calculate our estimation of the Monte Carlo expectation of the ratio of the /i-like 
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events to the e-like events. The numbers we obtain are within 10% of the values quoted 
by Kamiokande collaboration in Ref. 0. The differences could be due to the different set 
fluxes used and due to the simple approximation we made for the cross sections. 
In the presence of oscillations, the number of /i-like and e-like events are given by 



osc 



I [(t^,P,^cj + + C^eM^ + <i^,P-c-A (4c + 4c) dE (22) 



-''el osc 



blP,,a + <PlP,,a + ^i^P^^a + 4,Pp-ea\ e'dE (23) 



where P^/js are the probabilities for neutrino of flavor a to oscillate into flavor /3, derived 
in the last section. These oscillation probabilities are functions of the distance of travel 
(i, the mixing angles (j) and ip, the mass difference ^31 and the matter term A. These are 
calculated using the formulae in eqs. (|^), (^, (p!OD and (|T^). The probability of oscillation for 
anti-neutrinos P^^, in general, is different from P^p because of the different A dependence. 



From equations (p4 ) and (g^) we can calculate the ratio of //-like events to e-like events 
in the presence of oscillations to be 

rise = (^] (24) 



and the double ratio 



e / osc 



Ksc = -f^- (25) 

' MC 

If the atmospheric neutrino deficit is due to neutrino oscillations, then the double ratios 
Plsc given in eq. (^) should be within the range of the corresponding observed double ratios 
Rlijg, which are given in Table I. We searched for the values of the neutrino parameters 0, ijj 
and for which the predicted values of R^^^ satisfy the experimental constraints on the 
double ratios for all the five bins. The ranges of variation in the three parameters are 

1. < (f) < 50^. This is the range of (p allowed by the solar neutrino problem. For 
this range of (p, there exist values of ^21 and uj such that all the three solar neutrino 
experiments can be explained [|T^,^. 
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2. < ip < 90°. i/j is varied over its fully allowed range. 

3. 10~^ eV"^ < (5,31 < 10~^ eV^. The lower limit is given by the sub-GeV data and the 
upper limit is the largest value allowed by the two flavor analysis of the multi-GeV 
data by Kamiokande [Q. 

The results are plotted in Figures 2, 3 and 4. Figure 2 gives the projection of the allowed 
region on the 4> — ip plane, Figure 3 gives the projection on the — ^31 plane and Figure 4 
gives the projection on the ip ~ ^31 plane. The solid lines enclose the regions of parameter 
space whose predictions lie within the experimental range given by la uncertainties. The 
broken lines enclose regions whose predictions fall within range given by 1.6a uncertainties. 

As seen from Table I, the uncertainties in bin 5, which has {cos9) = 0.8, are quite large 
compared to the uncertainties in the other four bins. Morevoer, r^f,^ is greater than r^j^j 
through most of its range. Hence the double ratio -R^^^ > 1 for most of its range. The Monte 
Carlo expectation of the electron neutrino flux is less than that of the muon neutrino flux 
and the oscillation probabilities Pa/B are all less than 1. Using these facts, one can show from 
eqs. ( p!7D and ( p4D that, in general, r*^^ < r\.fQ or R^^^ < 1. Therefore, the region of overlap 
between R^f^g and R^g^ is very small. It is 0.9 — 1.0 for la uncertainties. It is possible that 
this small overlap is imposing a very strong constraint, leading to a situation where the bin 
with the largest uncertainty is essentially controlling the allowed values of the parameters. 
Because of this unsatisfactory situation, we redid the analysis ignoring the constraint from 
bin 5. We searched for regions of parameter space for which the values of i^^sc within 
the ranges of corresponding for only the flrst four bins. 

The results of the 4 bin analysis are plotted in Figures 5, 6 and 7. Figure 5 gives the 
projection of the allowed region on the (f) — ip plane. Figure 6 gives the projection on the 
<P~^3i plane and Figure 7 gives the projection on the 1^ — 631 plane. As before, the solid lines 
enclose the regions satisfying la vetoes and the broken lines enclose regions allowed by 1.6a 
vetoes. Comparing the corresponding flgures we flnd that the allowed values of parameters 
for the 4 bin flt are the almost identical to those from the 5 bin flt at 1.6a level. The allowed 
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regions at la level are somewhat larger compared to the 5 bin fit. This is not surprising 
because the bin 5 imposes a very strong constraint at Icr level. If this constraint is relaxed, 
then a somewhat larger region is allowed. The 4 bin fit shows that the 5th bin, which has 
the largest uncertainty, does not exercise undue influence on the selection of the parameter 
space. 

IV. DISCUSSION AND CONCLUSIONS 

The parameter space shown in Figures 1 to 4 , together with the allowed values for u and 
from our earlier work [|16|, provides a complete solution to the solar and the atmospheric 

neutrino problems in the framework of three flavor neutrino oscillations. The salient features 

of the results of the multi-GeV analysis are: 

• Most of the parameter space allowed by the multi-GeV data is a subset of the space 
allowed by the sub-GeV data. 

• The range of S31 allowed by la vetoes is extremely narrow. It is very close to the best 
fit value given by the two flavor analysis of Kamiokande. 

• The value of ip is always large {ip ^ 40°) and ip = 90° is allowed. 

• In the region allowed by la vetoes is always non-zero. = is allowed only at 1.6a 
vetoes. 

From the parametrization of , effective two level mixings can be obtained for the 
following choices of the angles: 

• Ue i^fj, for ip = 90°, 

• Uf, ^ i^r for ip = and 

• i^r ioT = 0. 
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Any solution to the atmospheric neutrino problem should suppress the muon neutrinos, 
enhance the electron neutrinos or do both. The v^, *^ channel, which suppresses electron 
neutrinos but leaves muon neutrinos untouched, cannot account for the atmospheric neutrino 
problem. Hence any solution of atmospheric neutrino problem should be away from the 
effective two flavor z/g ^ oscillations. The large value of the angle i\} is just a reflection 
of this fact. The allowed region includes the value ■?/' = 90°. Then the atmospheric neutrino 
problem is explained purely in terms of the two flavor oscillations between z/g v^^ with 
the relevant mass difference being ^31. In this case, the solar neutrino problem is solved by 
z/g — >• Vr oscillation which is determined by the mass difference ^21 (and the mixing angles lo 
and 0). 

How important are the matter effects in the analysis of the multi-GeV data? It can be 
seen from figures 2(d)-2(f) of Ref. that most of the expected multi-GeV events are caused 
by neutrinos with energies less than 10 GeV (over 80% for muon-like events and over 90% for 
electron-like events). The matter term, for a neutrino of energy 5 GeV, is about 2 x 10~'^e\^^. 
Since the initial range we considered for varied from 10~^ eV"^ to 0.1 eV^^, apriori one 
must include the matter effects in the expressions for the oscillation probabilities. However, 
the value of ^31 in the allowed region, especially for the la vetoes, where it is about 0.03 eV^^, 
is much larger than the matter term. Therefore, it is likely that the matter effects may not 
play an important role in determining the allowed parameter regions in the analysis of the 
multi-GeV data. To check this we reran our program with the matter term set equal to zero. 
With this change, the double ratios -Rq^c, defined in equation (p^), changes by about 10% in 
the first bin and by about 5% in the second bin. There is no discernible change in the other 
three bins. Since the errors in i?^^^ are about 30%, these small changes in i?^^^ do not lead to 
any apprecialble change in the allowed regions of the parameter space. However, the effect 
of matter terms may become discerible when more accurate data from Super Kamiokande 
become available. 

Since the earth matter effects seem to play no role in the determination of the parameter 
space, can one interpret the observed zenith angle dependence purely in terms of vacuum 
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oscillations? For an energy of 5 GeV, the mass square difference 10~^ eV"^ corresponds to 
an oscillation length of about 1200 Km. Thus bins 1 and 2 contain many oscillation lengths 
and the second and the third cosine terms in the vacuum oscillation probability, given in 
equation (|15[), average out to zero. For bins 4 and 5 the cosine terms are almost 1. In bin 3, 
these terms take some intermediate value. Therefore we have large suppression in the first 
two bins, almost no suppression in the last two bins and moderate suppression in the middle 
bin. 

Eventhough matter terms play no important role in the multi-GeV analysis, Ue I'n 
oscillations provide a better fit to data compared to z/^ ^ Ur oscillations. In ^ z/^ 
oscillations, the i/^ flux is suppressed whereas the flux is untouched. Thus the double 
ratios Rl^^ do become less than 1 but the large suppression observed in bins 1 and 2 is 
difficult to achieve via this channel. For z/g ^ oscillations, the z/^ flux is reduced and 
flux is increased. This occurs even when P^e = Pe/i, which is the case for vacuum oscillations, 
because the flux of z/^ is roughly twice the flux of z/g before oscillations. Hence the double 
ratios R^^^ are smaller for the case of z/g ^ z^^ oscillations compared to the case of z/^ z/^ 
oscillations and they fit the data better. 

We now make a brief comment on the LSND results on the search for z/^ — > z/g oscillations 



Tsj in the context of our analysis of atmospheric neutrinos. The LSND collaboration gives 



an oscillation probability P^g = (3.1_{q ± 0.5) x 10~ for muon anti-neutrinos in the energy 
range 20 — 60 MeV. In the framework described in section II, the oscillation probability 
relevant for the LSND experiment is the vacuum oscillation probability 



P',, = P'e = Sin' 20 sin^ ^ sin^ ( 1.27^ ) . (26) 



Note that both (p and ip have to be non-zero for P^g to be non-zero. In the region allowed 
by the Icr vetoes of the multi-GeV atmospheric neutrino data we have 

Minimum (sin^ 20sin^ ^ 0.04 for ~ 8°, ■j/' ~ 40° 
Maximum ^sin^ 20sin^ -0^ ~ 1 for ~ 40°, -0 ~ 90°. 

15 



Substituting these values and the oscillation probability obtained by LSND in eq. p6|), we 
obtain 

0.001 < sin^ ^1.27^^^ < 0.1. (27) 

For the LSND experiment the distance d = 30 meters. Taking the average energy to be 
(E) = 40 MeV, we obtain the range of ^31 to be 

0.03 eV^ < S31 < 0.3 eV^. (28) 

From the analysis of multi-GeV atmospheric neutrino data we have the upper limit on 
<^3i ^ 0.06 eV"^ (Figures 3 and 4). Hence there is a small region of overlap between the range 
of neutrino parameters required by the atmospheric neutrino data and the LSND data. This 



suggests that the standard three flavor analysis can accommodate all the data so far ||26 
and perhaps a fourth sterile neutrino is not needed. 

In conclusion, we have analyzed the atmospheric neutrino data of Kamiokande in the 
context of three flavor neutrino oscillations. We took into account both the zenith angle 
dependence of the multi-GeV data and the matter effects due to the propagation of the 
neutrinos through the earth. We obtained the regions in neutrino parameters which solve 
both the solar and the atmospheric neutrino problems. We found that the matter effects 
have negligible influence on atmopsheric neutrinos even in the multi-GeV range. The allowed 
regions of the parameter space with or without matter effects are almost identical. 



As we finished this work, we came across a preprint by Fogli et al ^7^. They have 
analyzed the sub-GeV data from various experiments and the zenith angle dependent multi- 
GeV data from Kamiokande. Although a detailed comparison is difficult, qualitatively our 
results agree with theirs. 
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TABLES 



Bin No. 


{cosO) 


[distance) in Km 




^obs 


Kbs 


1 


-0.8 


10,210 


3.0 




0.29l^:J? 


2 


-0.4 


5,137 


2.3 


l.UD_o 30 




3 


0.0 


<S32 


2.1 


-, 07+0.32 
-'-■^"-0.2:-! 


"■'JJ-o.ii 


4 


0.4 


34 


2.3 






5 


0.8 


6 


3.0 




1 O+0.6 



TABLE I. Zenith angle dependent data from Kamiokande [6] 
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FIGURES 

FIG. 1. Allowed region m. cj) — ip plane by the sub-GeV data (with 5zi > 10~^ eF^) at la 
(enclosed by solid lines) and at 1.6(T (enclosed by broken lines) 

FIG. 2. Allowed region vn (j) ~ ip plane by 5 bin analysis of multi-GeV data 
(with 10~^ eV"^ < 5z\ < 10~^ eV"^) at la (enclosed by solid lines) and at 1.6(7 (enclosed by 
broken lines) 

FIG. 3. Allowed region in 4> — ^3i plane by 5 bin analysis of multi-GeV data (with < < 90°) 
at la (enclosed by solid lines) and at l.Qa (enclosed by broken lines) 

FIG. 4. Allowed region in ^ — plane by 5 bin analysis of multi-GeV data (with < ^ < 50°) 
at la (enclosed by solid lines) and at l.Qa (enclosed by broken lines) 

FIG. 5. Allowed region in cf) — ij) plane by 4 bin analysis of multi-GeV data 
(with 10~^ eV"^ < S31 < 10 eF^) at la (enclosed by solid lines) and at 1.6a (enclosed by bro- 
ken lines) 

FIG. 6. Allowed region in cp — 631 plane by 4 bin analysis of multi-GeV data (with < ip < 90°) 
at la (enclosed by solid lines) and at 1.6a (enclosed by broken lines) 

FIG. 7. Allowed region intp — S31 plane by 4 bin analysis of multi-GeV data (with < cf) < 50°) 
at la (enclosed by solid lines) and at 1.6a (enclosed by broken lines) 
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